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Shapes and Sizes from Non-Identical-Particle Correlations 
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I review the prospects for measuring source ciiaracteristics from correlations other than those involving iden- 
tical pions. Correlations generated from Coulomb and strong interactions are shown to provide remarkable 
resolving power for determining three-dimensional information, in some cases accessing more detail than can 
be represented by Gaussian fits. 
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I. INTRODUCTION AND THEORY 

As is well-known to the participants of this conference, 
identical-particle correlations provide insight into the space- 
time structure of the asymptotic phase-space density iH. 
The experimentally measured two-particle correlation func- 
tion C(P,q) is related to the source function S (P,r) through a 
simple Fourier transform. 



3^(q)=Cp(q)-l = j (frSv{r)cQS.{q-r). 



(1) 



Here, P is the total momentum and q = (pi ~ p2)/2 is the 
relative momentum of the pair. The source function is defined 
by the outgoing phase space density. 



^p(r) 



Jd\id\2?>{ri n r)/(P/2,ri,f)/(P/2,r2,f) 
/jViJV2/(P/2,ri,f)/(P/2,r2,f) 



(2) 

where all quantities are calculated in the rest frame of the pair, 
where P = 0. As long as the time t is beyond any time for 
which particles are emitted, S is independent of t since both 
particles are moving with identical velocity, which is zero in 
this frame. The source function physically represents the nor- 
malized probability for two particles of the same velocity to be 
separated by a distance r in their asymptotic state. Although 
the source function only measures a distribution of relative 
distances, it still provides invaluable insight into the dynam- 
ics of a collision. For instance, long lived sources will lead to 
large separations along the direction of the pair's velocity. 

Determining S from the experimentally measured coiTela- 
tion function in Eq. ^ is simple as one needs only Fourier 
transform 2^ (q). Since the total momentum P is not involved 
in the convolution, I will suppress it in the notation for the 
remainder of the paper The ability of correlations of non- 
identical particles to reproduce source functions is less appre- 
ciated by the community, despite the fact that non-identical 
particle correlations offer similar resolvingpower to identical- 
particle correlations in many instances f£\. Non-identical- 
particle correlations are driven by the strong and Coulomb 
interactions. After making similar, though not identical, as- 
sumptions as what were used for Eq. ([TJ, the corresponding 
expression for the correlations are: 

^(q)=|«'V5(r)(|c^(q,r)|2-l). (3) 



Again, the label P is suppressed. If one views !^ (q) as a vec- 
tor, where each bin in q represents a component of the vec- 
tor, and 5(r) as a vector where each bin in r represents a 
component, determining S from R becomes a matrix inver- 
sion problem, where — 1 is the matrix with labels q and r. 
For the rest of the paper I will refer to the matrix as a kernel, 
3C(q,r) = |^|2-l. 

The goal of this talk is to understand the resolving power 
of the kernel. When viewed as a matrix, the kernel ?(;(q,r) 
is highly singular, i.e., it is not possible to uniquely determine 
every value of s{r) from even perfect knowledge of ?(,(q)- 
The information made accessible through the kernel depends 
on the unique properties of the wave function for a particular 
particle pair. For instance, the pA wave function is driven only 
by the strong interaction, which provides a large peak in (q) 
at small q. In contrast pK^ correlations are mainly driven 
by the Coulomb interaction, which tends to wash out the ef- 
fects of the strong interaction, but provides resolving power 
through the Coulomb interaction itself. Before discussing the 
kernels I first introduce the concept of expanding the prob- 
lem in spherical harmonics or Cartesian harmonics in the next 
section. In the subsequent sections, I review the general prin- 
ciples determining the resolving power of the kernel for both 
strong and Coulomb interactions. The resolving power for 
extracting non-Gaussian features a. pp source function is then 
presented as an example. 



II. THE RESOLVING POWER OF STRONG AND 
COULOMB INTERACTIONS 



Before launching into a detailed discussion of the resolving 
power of the kernel, 3C (q, r) = \^'\^ —1,1 project the coiTela- 
tion and source function with spherical and Cartesian harmon- 
ics 1 3]. The convenience of this decomposition derives from 
the fact that the specific projections of the correlations func- 
tion, ,„{q) for spherical harmonics or (,,,£- for Cartesian 
harmonics, are related to the same projection for the source 
function, 

^miq) = A-K. j r^dr St,n{r)'Ki{q,r), (4) 
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TABLE I: Cartesian harmonics for I <= 4. Other harmonics can 
be found by swapping indices on both sides of the equation, e.g., 
X y. For example, given ^210 = "?"v ^ iv/5, swapping y ^ z 
gives ^201 = — «j/5. 




^(4) 



(4) 
310 

„(4) _ 2 2_ 
■^220 " "x'^y 



(4) 
211 



:4-(6/7)n2 + (3/35) 
(l/7K2_(i/7)„2 + (l/35) 



For Cartesian harmonics, I = 
same for both expressions. 



£ + 4- The kernel is the 



1 



d cos 9„ 



/,r,cose^,-)| -1 Pi'(cose<,,-) 



(5) 

The projections reduce angular information F (£2) into a set of 
coefficients, F(„j or F^. 

The projections and inverses for spherical harmonics are 
defined as: 

Fe,„ = (47t)-i/2| daF{a)Ye^{a), (6) 

i,m 



The {2£ + 1 ) complex coefficients obey the symmetry F(,„ = 
(— so that the angular information for each £ is rep- 
resented by {2i+ 1) independent real numbers. 

For Cartesian harmonics, the projections are defined as: 



FiCl) 



(2^+1)!! fdD. 



£1 



4% 



F{a)A^{Q.), 



(7) 



^£J£y\£X(' 



Fjh^'^n/nl' 



The Cartesian harmonic function A^(£2) represent products of 

unit vectors, n/ny 'n^-, which then have lower values of £ pro- 
jected away. Table U shows the functions for £ < 4. Unlike 
the spherical harmonics, all these functions are real. For a 
given £, there are (^+ l)(^ + 2)/2 different combinations of 
{£x,£y,£z) that sum to £. However, not all coefficients are in- 
dependent. Since any part of the function from which one can 
factor n^ + n^ + n^ = 1 gives a function of lower projection £, 
the functions A-jt and the coefficients F-^ satisfy the "traceless- 
ness" constraint. 



Fe^+2/y/^+Fe,/^,+2/^+Fe,/^j^,+2 = 0. 



(8) 



With this constraint, all the coefficients for a given £ can be 
generated unambiguously from knowing the {2£+ 1) coeffi- 
cients for £x = or 1. Thus, both expansions, with spheri- 
cal harmonics or Cartesian harmonics, can be represented by 
{2£ + 1 ) real numbers for each £. 

The two choices of basis functions for the projections are 
equivalent, and it is straight forward to convert between one 



set of coefficients and the other. The main advantage of Carte- 
sian harmonics is that it is easier to identify the coefficient 
with the distortion of the shape as can be seen by the expan- 
sions in Eq. (|7]i. 



III. THE RESOLVING POWER OF 
IDENTICAL-PARTICLE STATISTICS 

The kernels for identical particles are particularly simple. 



i*C(q,r) = cos(q-r), 



(9) 



Inverting (q) to find S (r) involves a simple Fourier trans- 
form. Roughly speaking, the correlation structure of level q 
reveals information about components of the source function 
of scale R ^ 1/q. Shape characteristic are also easy to invert 
using the kernel for specific £. Since ji{qr) behaves as {qrY , 
information about the distortions of a specific £ require look- 
ing at a range of q for which qR > £. The fact that the kernel 
does not provide resolving power for small r is not usually a 
problem, since the source function Se{r) also tends to vanish 
for (r/R) «£. 



IV. THE RESOLVING POWER OF STRONG 
INTERACTIONS 

First, I consider the information available by analyzing the 
kernel for £ ^0. For q >> 1/R, where R is a characteristic 
scale of the source, |(|)(^,r)p differs significantly from unity 
only at small r, and the kernel can be approximated by a delta 
function. The correlation function is then only sensitive to the 
probability that the two particles have r = 0, with a resolving 
power determined by the overall density of states which is 
given by the scattering phase shifts. 



Adn /dq 



Adn/dq 



dq 



(10) 
(11) 



The sensitivity in q is determined by the form of d?> / dq, rather 
than by details of the shape. 

For ^ < 1 /R, where /? is a characteristic scale of the source, 
^=o{q) becomes sensitive to more details of Se=Q{r), not just 
the value at r = 0. At that point, the resolving power depends 
on the unique shape of the wave function, which is also given 
by the phase shifts. For ^ << 1 /R, the kernel is determined 
solely by the £ = Q phase shifts. As long as the range of the 
interaction is short compared to the source size (true for heavy 
ion collisions), one can then approximate the wave function by 
a distorted plane wave. If only the £ — pieces are kept. 



«c(q,r) 



sin(q'r + 5^=0) sin<7r 
qr 

3C(q^0,r): 



1 



qr 

—2ar + a^ 



(12) 
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where a is the scattering length, d(=Q = —qa. Thus, one can 
crudely state that the %ji=Q{q 0) provides information about 
the and moments of the source function, while 

%i=(i{q » 1/^) provides information about Si=(j{r 0). 
Looking at d%£=Q jdq^ at ^ will provide insight into other 
higher moments such as ( 1 ) and (r) . For higher moments the 
constant of proportionality depends on the l—\ phase shifts 
and on the next higher expansion of the £ = phase shift, 
which is determined by the effective range. 

Gaussian parameterizations of Si={){r) involve two param- 
eters, a size Rq, and a "coherence" parameter A, which repre- 
sents the normalization of S which can differ from unity since 
some of the particles might originate from effectively infinite 
distances due to long-lived decays. Any two pieces of infor- 
mation about Si=o{r), such as the value at r = and a moment, 
are sufficient to determine the two Gaussian parameters. As 
explained in the previous paragraph, the kernel can, in princi- 
ple, provide more information. However, the ability to extract 
such non-Gaussian details are constrained by the strength and 
form of the scattering phase shifts and by statistical and sys- 
tematic errors in the data. 

The kernel for ^ > allows one to determine properties of 
the source shape. Even though the scattering phase shifts are 
usually important only for a few small values of I, %i{q, r) has 
strength for all £. This sensitivity can be explained by pointing 
out that for large qr, which is the only region one can analyze 
higher I, one can view the problem classically. The kernel 
is then largely determined by the shadowing. In the large qr 
limit the kernel can be determined solely from the geometry 
of classical trajectories. 



(q,r) 1.^00 = - 



1 f c . d(3 
;3 -08(-il„-) + ^ 



(13) 



Here, the first term represents the shadowing due to particles 
scattering away from the direction q-f — ~\, while the second 
term represents the addition due to scattered particles. Since 
the delta function has contributions to all the kernel from 
strong interactions provides resolving power at all ^ for larger 
r. In practice, the power falls off at large t since quantum 
considerations destroy for qr < I. 

The t = Q kernel shown in Eq. (fTOl i and the ^ > kernel 
shown in Eq. ( flJl l tend to be of similar strength for most cir- 
cumstances. Thus, for most interactions, it can be stated that 
if one has sufficient statistics to measure the size, one also has 
the ability to determine shape. 



V. THE RESOLVING POWER OF COULOMB 
INTERACTIONS 



to the kernel. For most experiments the resolution is of the or- 
der of 1 MeV/c and only those pairs where at least one mem- 
ber originates from a weak decay are uncorrected. Thus, the 
intercept of the correlation function at ^ = gives 



C(<?^0) = (l-/i,.)2 



(14) 



where /i.i. is the fraction of particles from long-lived decays. 

For large q, one can use the classical limit of the kernel 
given by calculating classical Coulomb trajectories. 



|(j)(^,r,cose^r)P ^ 

1 +C0s9^r — -^^ 



(15) 



=©(1 +cose,,,.-2x). 



^(1 +C0s9f/r--«)^ 

Here, x is the ratio of the Coulomb and final-state kinetic en- 
ergy, 



q^r 



(16) 



Since = 1/137 is small, x is small unless q is quite small. 
This allows the kernel to be approximated as 

aC(q,r)|,^o « -^5(l+cos(e,,), (17) 

In addition to x being small, this approximation also relies on 
the classical approximation, qr >> £. 

The Coulomb interaction thus provides excellent resolving 
power for all £, provided that x is not too small for the rel- 
evant region qR > I, which can equivalently be stated that 
/je^R should not be too small, or the source size should not be 
too much smaller than the Bohr radius, which is 390 fm for tot 
but only 58 fm for pp. Thus, one expects reasonable resolving 
power for pK or heavier pairs, whereas pairs involving pions 
offer weak resolving power through the Coulomb force. The 
Coulomb interaction becomes especially useful if light nuclei 
can be used, as in addition to the increase in the reduced mass 
/J, the increased factor of charges, Z\Z2e^, bolsters the 
resolving power 

For large r the Coulomb kernel falls as (^^r)^' for large q. 
Thus, the tail of the correlation function provides a measure of 
the moment (1/r) for the source function. By analyzing the 
correlation function for its variance from the behavior 
for decreasing q, one can, in principle, ascertain other facets 
of the source function. 



Coulomb interactions also provide leverage for determin- 
ing source sizes. The limitations of the Coulomb interaction 
derives from the inherently weak couplings which make corre- 
lations relatively small except at small q. At small q, 0, 
unless is so large as to be outside the range of the Coulomb 
interaction, e^/r << q^ /2^. Thus, all pairs are emitted within 
a radius, r < /je^/q^, where qr is the momentum resolution, 
contribute zero to the correlation function, or equivalently — 1 



VI. RESOLVING POWER IN THE REAL WORLD 

Imaging has become a popular term for describing the in- 
version of ^{q) to determine S{q)- Particular attention has 
been given to the ^ = cornponents by Danielewicz and 
Brown by fitting to splines flsla], where non-Gaussian as- 
pects of source functions have been observed. The spline- 
fitting routines are themselves a fit to a constrained functional 
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form, and should not be mistaken for a fully flexible account- 
ing for all possible forms for the source function. Our goal 
here is to gain a better feel for the sensitivity to non-Gaussian 
features of the source one might attain in real-world analyses. 

First, I will consider the simple case of £ = sources. I pick 
the form 



Se=Q{r) = 'k{l -Xfrac)- 



frac " 



(18) 

Here, the normalization constant ensures that each source 
would integrate to unity when A, = 1. The four parameters 
in this fit are: 

• A,, the "coherence factor". 

• Rg, a Gaussian source size. 

• X, an exponential size describing the source function at 
large r. 

• ^frac, the fraction of the source described by the expo- 
nential. 

The form of the argument of the second term was chosen to 
force the curvature at the origin to match that of the Gaussian. 
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FIG. 1: A one-dimensional pp correlation function for the param- 
eters shown in Eq. dlSt is displayed in the upper panel, while 
changes corresponding to the labeled changes of the four parame- 
ters are shown in the lower panel. For most experiments, only three 
of the four parameters would be uniquely determined. 

The upper panel of Fig. [T] shows a correlation function 
calculated for the form above with the parameters: A, = 0.6, 
Rq ~ 5 fm, X = \Q fm, Xf^c — 0.5. In the lower panel, the 
change of the correlation function is illustrated for various 
changes in the four parameters above. Visually, it is apparent 
that any one of the four changes can be pretty well approxi- 
mated by making linear combinations of the other three. For 



instance, the effects of increasing the exponential length X by 
2 fm can be approximately reproduced by a combination of re- 
ducing A combined with small changes to the remaining two 
parameters. 

The ability of an experimental analysis to determine all four 
parameters in Eq. ( fTSl l depends on details of the particular in- 
teraction and the experimental accuracy and resolution. To 
more rigorously determine the degree to which parameters 
could be experimentally constrained, theoretical correlation 
functions were binned, and random errors were then added. 
For errors of the order of what might be expected at RHIC, 
it did not seem possible to uniquely determine all four pa- 
rameters. However, extracting three parameters was usually 
quite robust. Since a Gaussian parameterization uses only two 
parameters, our exercise demonstrates that non-Gaussian fea- 
tures can be extracted, but that determining both the weight 
Xfiac and the exponential scale X will probably be outside the 
resolution of the experiments. However, if one feels confident 
in the scale X, e.g., one might believe it comes principally 
from (0 decays, our investigations suggest that it is quite pos- 
sible to confidently determine the relative weight Xfiac. Most 
of the pairs that were investigated appeared to accommodate 
three-parameter analyses. 
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FIG. 2: Correlations for two relative momenta are shown as a func- 
tion of the angle with respect to the z axis for a Gaussian source of 
radius R-^ = Ry = 4, R, = 8 fm. The angular variations for q = 55 
MeV/c illustrate significant resolving power for experimentally de- 
termining shapes. 
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As shown in the previous section, insight into the three- 
dimensional shape of sources can come from all types of in- 
teractions. This is illustrated by the plotting of several angular 
correlations for a non-spherical source in Fig. |2] The Gaus- 
sian source has dimensions of 4 fm in the x and y directions 
and 8 fm in the z direction. For non-identical particles, the 
center of the Gaussian was moved by 4 fm to simulate one 
species being emitted ahead or behind of the other ill S, |^. 
The correlation function is plotted at two values of as a func- 
tion of the angle with respect to the z axis. The fact that the 
values are not constant, and that their variation with angle is 
similar to their average strength, illustrates the strong angular 
resolving power provided by all the pairs. In all cases, their 
is little sensitivity to shape at small q. This can be under- 
stood from considering the uncertainty principle that says that 
if qR < 1, there is an inherent vagueness as to the direction 
of the relative momentum. However, at a relative momen- 
tum of 55 MeV/c, the angular sensitivity is significant for all 
cases. The pairs involving two baryons tended to be espe- 
cially strong. To illustrate the same information with angular 
projections. Fig. [3] shows projections in terms of Cartesian 
harmonics for all ^ < 3. 
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FIG. 3: Angular projections of correlations as a function of q for a 
Gaussian source of dimensions =Ry = 4 fm, /?z = 8 fm. For non- 
identical particles, the center of the Gaussian was also moved by 4 
fm in the z direction. Projections are for Cartesian harmonics, with 
£ = (lower panels), £x = iy = 0,(.z = li2,3, subsequently higher 
panels. 

As a point of comparison, I present identical-pion correla- 
tions for the same source, (Rx = Ry = 4 fm, .Rj = 8 fm), but 
without the offset of the center due to the fact that there can 
be no offset with identical particles. As can be seen in Fig. |4] 
the correlations are much larger than those seen for any of the 



examples from Fig. |3] The stronger resolving power derives 
from the simple fact that all the correlations, including those 
for £ = Q, are several times larger For identical particles, there 
are no odd-numbered harmonics. However, as shown in Fig. 
m it should be feasible to measure harmonics to rather large i. 
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FIG. 4: Angular projections of identical-pion correlations as a func- 
tion of q for a Gaussian source of dimensions = Ry = 4 fm, Rz = S 
fm. Since the correlations are stronger for this case, larger values of 
£ can be explored. 

In addition to the calculations for Gaussians shown here, I 
have evaluated angular projections for correlations from blast 
wave models. From our experience thus far I can make the 
following rough statements regarding the prospects for dis- 
cerning shapes. 

• For identical pion or identical kaon correlations, many 
values of £ are accessible. 

• For pK or pK correlations: 

- £ = I terms are easy 

- £ = 2 require 1 % or better accuracy 

- ^ = 3 are probably prohibitively small 

• For baryon-baryon correlations 

- £ = 1,2 are easy 

- £ = 3 requires 1 % or better accuracy. 



VII. SUMMARY AND PROSPECTS 

Identical-pion correlations arguably represent the most dis- 
criminating observable at RHIC. A large fraction of the space 
of possible equations of state have been excluded by the- 
ory/experiment comparison, including strong first-order equa- 
tions of state, and on the opposite extreme, hard pion-gas-like 
equations of state. Given the paramount importance of these 
conclusions, it would be enormously important to verify these 
statements with independent analyses of space-time parame- 
ters using other classes of correlations. As a first goal, exper- 
imental analyses should work on alternative measurements of 
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^out/^side/^iong- Additionally, the odd harmonics access new 
features of the emission geometry. These alternative mea- 
surements should be viewed as being more complementary 
than redundant, as the theoretical basis for Coulomb-induced 
and strong-interaction based correlations is somewhat differ- 
ent from that used for identical-particle interference. 

Exploiting other classes of shape analyses is inherently dif- 
ficult for two reasons. First, the strength of the correlations 
tends to be small. Whereas identical pion correlations are 
measured in the tens of percent, other classes of correlation 
tend to be a few percent once they are in the region <7 > 25 
MeV/c, where one can access shape information. Recent 
data sets from RHIC have amassed sufficient statistics to an- 
alyze correlation functions at the sub-one-percent level. The 
PHENIX experiment, with its excellent particle-identification 
over a wide range of momenta, is especially well poised to 
exploit such correlations. STAR has even better statistics 
for some cases, but the delay in finishing the STAR time-of- 
flight wall has made it difficult to access some correlations at 
high statistics. For instance, since correlations are analyzed at 
small relative velocity (not small relative momentum) p% cor- 
relations require comparing a p, = 100 MeV/c pion to a 700 



MeV/c proton to match velocities. Both particles are at the 
edge of STAR'S acceptance for particle identification. How- 
ever, analyses will ultimately be constrained by systematic 
experimental errors, and from competing physical processes, 
such as flow or jets. All such effects must be handled more 
carefully when the desired correlation is of order one percent. 

The second challenge for analyzing these classes of cor- 
relation derives from the difficulty in disentangling the more 
complicated and opaque kernels. However, as was shown 
here, this is by no means insurmountable. Recently devel- 
oped fitting algorithms, imaging techniques, and angular de- 
compositions, have now been compiled into usable libraries 
with the CorAL (Correlations Analysis Library) project. An 
alpha version of the library is already available [10], and in- 
terested parties are encouraged to contact any of the authors, 
for assistance in using the codes. 
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I. FORMAT 

The very first thing to do is to make sure you have RevTeX 
4 in the subdirectory where you will edit your contribution. In 
case you do not have it yet, please download the package form 
the internet and install it in our subdirectory. 

The text of your manuscript should be produced using 
two-column format, within the dimensions shown on these 
pages. The ETgX document class uses the maximum stipu- 
lated length apart from the following two exceptions (i) I^TgX 
does not begin a new section directly at the bottom of a page, 
but transfers the heading to the top of the next page; (ii) lATgX 
never (well, hardly ever) exceeds the length of the text area 
in order to complete a section of text or a paragraph. Some 
examples of reference citation follows, with the references at 
the end, as given by Ref. [? ? ] and in [? ]-[? ]. 



A. Spacing 

We normally recommend the use of 1 .0 (single) line spac- 
ing. However, when typing complicated mathematical text 
I^Tj5C automatically increases the space between text lines in 
order to prevent sub- and superscript fonts overlapping one 
another and making your printed matter illegible. 



A. Examples 

Tables should be presented in the form shown in Table ?? 
(when an intermediate line through all the columns is neces- 
sary) or in Table ??. They also illustrate cases with differ- 
ent number of columns. Their layout should be consistent 
throughout the text, so please use these models adapted to 
your contents. 



TABLE I: Table format with an intermediate line and 3 columns 





No Resonances 


LUND Resonances 


Comaparing Data to Model - 2 scenarios 


lxL-i|/<^ 


2.1 


2.2 




4.6± 0.9 


3.1± 1.3 


Ato 


3.4 ± 0.7 


1.6± 1.0 


a 


0.027 


0.014 


P 


0.042 


0.023 



The above example illustrates the case where mainly num- 
bers are included. In the example that follows, we show equa- 
tions in the table, which now contains less columns than the 
one above. 



B. Fonts 

These instructions have been produced using 11 point 
Times fonts, normally used by RevTeX 4 document class. 
Don't use too small fonts to guarantee legibility. 



II. TABLES 

Tables should be made with lATgX; illustrations should be 
originals or sharp prints. All the captions should be clearly 
legible. They should be arranged throughout the text and 
preferably be included on the same page as they are first dis- 
cussed. They should have a self-contained caption and be 
aligned in the center of the text. 



TABLE II: Table format with equations and 2 columns 



Source p(r,f) 


C{h,k2)-l 


A: 




B:5(r-i;) e-'Vz-t' 


g-(Aa))V [s\ri{qR)l{qR)f/2 


C:&(R-r) e-'^l^^- 


(9/2) e-(^")''^'{(<7i?)}-4 
{[coi{qR)-^m{qR)l{qR)]f 


D: e-'l'^ 0(3x2 2) 


[sin(AcoV3T)/AcoV3T] 
{l+q^R^)-'/2 


E: 

®(Rt-r)e-''l^''@{t) 


9|/|V(8a'''), 

I=-iy/Ti[{l^lJZ+)W{z+)- 

(\-tiz-)W{z-)]-2^ 



The captions to the tables should be included at the top, as 
shown in both Table ?? and Table ??. 



